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ABSTRACT 

We present new, high resolution, infrared spectra of the T dwarf Gliese 229B in the J, H, and 
K bandpasses. We analyze each of these as well as previously published spectra to determine 
its metallicity and the abundances of NH3 and CO in terms of the surface gravity of Gl 229B, 
which remains poorly constrained. The metallicity increases with increasing gravity and is below 
the solar value unless Gl 229B is a high-gravity brown dwarf with log<?(cgs) w 5.5. The NH3 
abundance is determined from both the H and the K band spectra which probe two different 
levels in the atmosphere. We find that the abundance from the K band data is well below 
that expected from chemical equilibrium, which we interpret as strong evidence for dynamical 
transport of NH3 in the atmosphere. This is consistent with the previous detection of CO and 
provides additional constraints on the dynamics of the atmosphere of this T dwarf. 

Subject headings: stars: individual (Gl 229B) — stars: brown dwarfs — stars: atmospheres — stars: 
abundances 
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1. Introduction 

Gliese 229B is not only the first brown dwarf 
recognized as genuine (Nakajima et al. 1995; Op- 
pcnheimcr ct al. 1995), but it is also the brightest 
and best-studied T dwarf known. With an effec- 
tive temperature of T e g ~ 950 K, it lies squarely 
between the latest L dwarfs (T e g ~ 1500 K, Kirk- 
patrick et al. (1998)) and the giant planets of the 
solar system (T e g ~ 100 K). Indeed, its near in- 
frared spectrum shows the strong H 2 absorption 
bands characteristic of very- low mass stars and the 
strong CH4 bands seen in the spectra of Jupiter, 
Saturn and Titan. The transitional nature of the 
spectrum of Gl 229B is remarkable and hints at 
the spectral appearance of extrasolar giant plan- 
ets which have effective temperatures in the range 
200 - 1600 K (Guillot 1999). 

A wealth of data on Gl 229B has accumulated 
since its discovery five years ago. Broad band 
photometry from R through N and an accurate 
parallax (Matthews et al. 1996; Golimowski et al. 
1998; Leggett ct al. 1999; Perryman ct al. 1997) 
allow an accurate determination of its bolometric 
luminosity. Spectroscopic observations (Oppen- 
hcimer et al. 1998; Geballe et al. 1996; Schultz 
ct al. 1998) covering the range from 0.8 to 5.0 ^m 
have revealed a very rapidly declining flux short- 
ward of 1 /im, the unmistakable presence of CH4, 
H 2 0, and Cs, and demonstrated the absence of the 
CrH, FcH, VO and TiO features characteristic of 
late M and early L dwarfs. Finally, Noll, Geballe & 
Marley (1997) and Oppenheimer et al. (1998) have 
detected CO with an abundance well above the 
value predicted by chemical equilibrium, a phe- 
nomenon also seen in the atmosphere of Jupiter. 

Model spectra for Gl 229B (Marley et al. 1996; 
Allard et al. 1996; Tsuji et al. 1996b) reproduce 
the overall energy distribution fairly well and all 
agree that 1) T e g ~ 950 K, 2) compared to gaseous 
molecular opacity, the dust opacity is small if not 
negligible in the infrared, 3) the gravity of Gl 229B 
is poorly constrained at present. The rapid de- 
cline of the flux at wavelengths shortward of 1 /zm 
is interpreted as caused by an absorbing haze of 
complex hydrocarbons (Griffith, Yelle & Marley 
1998) or alternatively by the pressure-broadened 
red wing of the K I resonance doublet at 0.77 ^m 
(Tsuji, Ohnaka & Aoki 1999; Burrows, Marley & 
Sharp 1999). 



In this paper, we present new high-resolution 
spectra in the J, H, and K bands. With the inclu- 
sion of the "red" spectrum of Oppenheimer et al. 
(1998), we analyze each part of the spectrum sepa- 
rately to obtain independent measures of the H 2 
abundance of Gl 229B - broadly interpreted as the 
metallicity index - to detect for the first time the 
presence of NH 3 in its spectrum, and to estimate 
the NH 3 abundance at two different depths in the 
atmosphere. Our results are expressed in terms 
of the surface gravity which cannot be determined 
from the data presented here. Nevertheless, we 
identify a reduced set of acceptable combinations 
of T e g and gravity, using the observed bolometric 
luminosity of Gl 229B (Leggett et al. 1999). 

The observations and the near infrared spec- 
tra are presented in §2. Section 3 shows how an 
accurate parallax, a well-sampled spectral energy 
distribution and evolutionary models greatly re- 
duce the possible range of combinations of T e g and 
gravity without having to resort to spectrum fit- 
ting. The synthetic spectrum calculation and our 
method of analysis are described in §4. The results 
concerning several molecules of interest which are 
at least potentially detectable are presented in §5, 
followed by a discussion in §6. Finally, a summary 
of the results and directions for future study are 
given in §7. 

2. Spectroscopic observations 

Spectra of Gl 229B in selected narrow inter- 
vals in the J, H, and K windows were obtained 
at the 3.8 m United Kingdom Infrared Telescope 
(UKIRT) in 1998 January, using the facility spec- 
trometer CGS4 (Mountain et al. 1990) and its 150 
1/mm grating. Details of the observations are pro- 
vided in Table 1. These are among the highest 
resolution spectra obtained of any T dwarf. 

The spectra were obtained in the standard 
stare/nod mode with the 1.2" wide slit of the spec- 
trometer oriented at a position angle of 45° , nearly 
perpendicular to the line connecting Gl 229A and 
Gl 229B. The southward-going diffraction spike of 
Gl 229A together with scattered light from that 
star, which is 10 magnitudes brighter than Gl 
229B, contaminated the array rows near and to the 
southwest of those containing the spectrum of Gl 
229B. The contamination on the Gl 229B rows was 
determined by interpolation and was subtracted; 
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typically it was comparable or somewhat smaller 
than the signal from Gl 229B. In order to remove 
telluric absorption features, spectra of the AOV 
star BS 1849 were measured just prior to Gl 229B. 
In all cases the match in airmasses was better 
than five percent and hence in the ratioed spectra 
residual telluric features are small compared to the 
noise level. Wavelength calibration was achieved 
by observations of arc lamps and is in all cases 
better than one part in 10 4 (2a). 

The spectra shown in this paper have been 
slightly smoothed, so that the resolving powers are 
lower than those in Table 1 by approximately 25 
percent. They also have been rebinned to facilitate 
coadding like spectra and joining adjacent spec- 
tral regions. The error bars can be judged by the 
point-to-point variations in featureless portions of 
the spectra, the signal-to-noise ratios at the con- 
tinuum peaks are approximately 40 in the K band, 
25 in the H band, and 30 in the J band. The flux 
calibration of each spectrum is approximate as no 
attempt was made to match the photometry of Gl 
229B. 

While we identify the spectra by their cor- 
responding standard photometric infrared band- 
pass, their wavelength coverages are much nar- 
rower than the JHK filters and typically corre- 
sponds to the peak flux of Gl 229B in each band- 
pass. 

In the J band spectrum (Fig. 1), nearly all fea- 
tures are caused by H 2 0. The short wavelength 
end of the spectrum shows the red side of a CH4 
band, which is responsible for the features seen 
shortward of ~ 1.215 /im. Two lines of neutral 
potassium are easily detected near 1.25 /im. No 
other lines of alkali metals fall within the wave- 
length coverage of our observations. The H band 
spectrum (Fig. 2) is relatively rich in molecu- 
lar opacity sources. All the features seen in the 
spectrum are either due to H 2 (A & 1.59 /im) 
or part of a very strong CH4 absorption band 
(A ^ 1.6 ^m). Features seen between 1.59 and 
1.6 /im cannot presently be ascribed with certainty 
and are due to either H2O or CH4. While the opac- 
ities of NH 3 and H 2 S are not negligible in this part 
of the spectrum, neither molecule forms distinctive 
spectral features. Their presence cannot be di- 
rectly ascertained from these data, mainly because 
their opacity is weaker than that of H 2 and CH 4 
and because of significant pressure broadening (see 



§5.3). The K band flux emerges in an opacity win- 
dow between a strong H 2 band and a strong CH 4 
band (Fig. 3). Spectral features are caused by 
H 2 at shorter wavelengths (A < 2.11 /im) and by 
CH4 at longer wavelengths (A > 2.105 /im). Sev- 
eral distinctive features of NH 3 are expected at 
the blue end of this spectrum and models predict 
a single absorption feature of H 2 S at 2.1084 /im. 

All features seen in Figures 1-3 are unresolved 
blends of numerous molecular transitions. A spec- 
tral resolution at least 10 times higher would be 
required to resolve the intrinsic structure of the 
spectrum of Gl 229B. 

3. Effective temperature and gravity 

While synthetic spectra have been fairly suc- 
cessful at reproducing the unusual spectrum of 
Gl 229B (Marley et al. 1996; Allard et al. 1996; 
Tsuji et al. 1996b), the entire spectral energy 
distribution has not yet been modeled satisfac- 
torily. Limitations in the opacity databases are 
partly responsible for the remaining discrepan- 
cies between synthetic spectra and the data (see 
§4.2). These shortcomings have impeded the de- 
termination of T e g and of the gravity g in par- 
ticular. On the other hand, the bolometric lumi- 
nosity of Gl 229B is now well determined. Com- 
bining spectroscopic and photometric data from 
0.82 to 10 /im with the parallax, Matthews et 
al. (1996) found L = 6.4 x lO" 6 !,©. With new 
JHKL' photometry, Leggett et al. (1999) found 
L = 6.6 ± 0.6 x 10~ 6 L Q . A recalibration using 
the HST photometry of Golimowski et al. (1998) 
gives L = 6.2 ± 0.55 x 1O~ 6 £ . Evolutionary 
models (Burrows et al. 1997) allow us to find a 
family of (Tcffjff) models with a given Lboi- Fig- 
ure 4 shows the cooling tracks of solar metallicity 
brown dwarfs in terms of the surface parameters 
T c ff and g. Models with the bolometric luminosity 
of Gl 229B fall within the band running through 
the center of the figure. Using a very conservative 
lower limit for the age of Gl 229A of 0.2 Gyr (Naka- 
jima et al. 1995), we find T cS = 950 ± 80 K. On 
the other hand, the gravity remains poorly deter- 
mined with \ogg (cm/s 2 ) = 5 ± 0.5, corresponding 
to a mass range of 0.015 - 0.07 M . While the up- 
per range is very close to the lower main sequence 
mass limit, Gl 229B's status as a brown dwarf is 
secure. A star at the edge of the main sequence 
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would be much hotter with T e g- ~ 1800 K; well 
outside of Fig. 4. 

In the remainder of this paper, the discussion 
focuses on three atmosphere models which span 
the range of allowed solutions (Table 2 and Fig. 
4): (T cff (K),log.g(cgs)) - (870, 4.5), (940, 5.0) 
and (1030, 5.5), which we label models A, B, and 
C, respectively. These constraints on T c g and g 
from cooling sequences are quite firm. We find the 
same result, within the error bar on -Lboii from sev- 
eral cooling sequences which predate Burrows et 
al. (1997). The latter were computed with differ- 
ent input physics such as the equation of state and 
surface boundary conditions derived from grey and 
non-grey atmosphere models using several opacity 
tabulations. 

In section 5, we show that these three mod- 
els can fit the spectra only if they have differ- 
ent mctallicities, ranging from [M/H]= —0.1 to 
—0.5. The evolution of brown dwarfs is sensitive 
to the metallicity through the atmospheric opacity 
which controls the rate of cooling. We find that 
for the range of interest here, the effect of a re- 
duced metallicity on our determination of T e ff, <?, 
and the cooling age is smaller but comparable to 
that of the uncertainty on the value of L^ \. We 
choose to ignore it for simplicity. 

4. Method of analysis 

4.1. Model atmospheres and spectra 

Our analysis is based on the atmosphere mod- 
els of brown dwarfs and extrasolar giant planets 
described in Burrows et al. (1997). Briefly, the at- 
mospheres are in radiative/convective equilibrium 
and the equation of radiative transfer is solved 
with the k-coefficient method. The chemical equi- 
librium is treated as in Burrows et al. (1997). Gas 
phase opacities include Rayleigh scattering, the 
collision-induced opacity of H 2 and the molecular 
opacities of H 2 0, CH 4 , NH 3 , H 2 S, PH 3 , and CO, 
as well as the continuum opacities of H~ and . 
The molecular line opacity database is described 
in more detail in §4.2. Atomic line opacity is not 
included. Because of the relatively large gravity 
of Gl 229B, pressure broadening of the molecular 
lines plays an important role in determining the 
(T, P) profile of the atmosphere and in shaping 
the spectrum. The line-by-line broadening theory 
we use is described in Burrows et al. (1997). The 



strong continuum opacity source responsible for 
the rapid decrease of the flux of Gl 229B shortward 
of 1.1 /im is included following the haze model of 
Griffith et al. (1998). Details of the haze model 
and of our fitting procedure are given in §5.1. The 
(T, P) structures of these atmosphere models are 
shown in Fig. 5. The profiles intersect each other 
because both T e g and the gravity vary between 
the models. The inflexion point at logT ~ 3.25 
signals the top of the convection zone. 

Using the same monochromatic opacities used 
to compute the k-coefficients, high-resolution syn- 
thetic spectra are generated from the atmospheric 
structures by solving the radiative transfer equa- 
tion with the Feautrier method on a frequency 
grid with Av — 0.1 cm -1 . Spectra with resolution 
lower than <; 1 cm -1 can then be generated 
for comparison with data. 

An unusual aspect of T dwarf atmospheres is 
the great variation of the opacity with wavelength. 
These atmospheres are strongly non-grey and the 
near infrared spectrum is sculpted by strong ab- 
sorption bands of H 2 and CH4. Most of the flux 
emerges in a small number of relatively transpar- 
ent opacity windows. The concept of photosphere 
becomes rather useless since the level at which the 
spectrum is formed depends strongly on the wave- 
length. Figure 6 shows the depth of the "pho- 
tosphere" [t v = 2/3) in both temperature and 
pressure as a function of wavelength for model 
B (X^ = 940 K, log.g = 5). In the Z, J, and 
H bands, and, to a lesser extent, in the K and 
M bands, the atmosphere is very transparent and 
can be probed to great depths. For A <; 25 fim, the 
spectral energy distribution approaches a Planck 
function with T ~ 500 K. 

Figure 6 reveals that spectroscopy between 0.8 
and 12 jim. can probe the atmosphere from T ~ 
500 K down to a depth where T ~ 1600 K, cor- 
responding to a range of 6 pressure scale heights! 
This provides an exceptional opportunity to study 
the physics of the atmosphere of a brown dwarf 
over an extended vertical range. The top of the 
convection zone for model B (located at T = 
1860 K) is below the "photosphere" at all wave- 
lengths and is not directly observable, however. 
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4.2. Limitations to this study 

Given the range of acceptable values of T e ff and 
g (§3), we can determine the metallicity of Gl 229B 
and the abundance of several key molecules by fit- 
ting synthetic spectra to the observations, for each 
of the three models. The precision of our results 
is determined by the reliability of the models, the 
noise level in the data and, most significantly, by 
the limitations of the molecular line lists used to 
compute their opacities. The latter point requires 
a detailed discussion. 

The opacities of CH4 and NH 3 are computed 
from line lists obtained by combining the HITRAN 
(Rothman et al. 1998) and GEISA (Husson et al. 
1994) databases, which are complemented with re- 
cent laboratory measurements and theoretical cal- 
culations. Further details are provided in Burrows 
et al. (1997). The resulting line lists for these two 
molecules are very nearly complete for T < 300 K, 
and their degree of completeness decreases rapidly 
at higher temperatures where absorption from ex- 
cited level become important. Furthermore, the 
line list of CH4 is limited to A > 1.584 /im. We 
extend the CH4 opacity to shorter wavelengths 
with the laboratory measurements of Strong et 
al. (1993) which provide the absorption coefficient 
averaged over intervals of 5cm _1 between 1 and 
5/xm at T = 300 K. We use Strong et al. (1993) 
opacities for 1 < A < 1.82 /im and the line list for 
A > 1.82 /im. This puts the transition from one 
tabulation to the other in a strong H2O absorp- 
tion band and obliterates any discontinuity in CH4 
opacity at the transition. For A < 1 /im, the tabu- 
lation of Karkoschka (1994) gives the absorption 
coefficient of CH4 determined from spectroscopic 
observations of the giant planets at 0.0004 /im in- 
tervals. Because of the low temperatures found in 
the atmospheres of giant planets, the Karkoschka 
CH4 opacities are appropriate for T ^ 200 K. To 
our knowledge, this compilation of NH 3 and CH4 
opacity is the most complete presently available. 

As can be seen in Figure 6, the temperature in 
the atmosphere of Gl 229B is everywhere greater 
than 300 K. For CH4, we compute temperature- 
dependent line opacity (which is incomplete above 
300 K) from the population of excited levels deter- 
mined by the Boltzmann formula for A > 1.82 /tm, 
and use temperature independent opacity at 
shorter wavelengths. While the synthetic spec- 



tra computed reproduce the fundamental band of 
CH4 very well (centered at A = 3.4 /tm), the match 
with the 1.6 and 2.3 /im bands is rather poor. Even 
though CH4 is a very prominent molecule in the 
spectrum of Gl 229B, the current knowledge of its 
opacity is not adequate for a quantitative analy- 
sis of its spectral signature. For this reason, we 
have essentially ignored the regions in our spectra 
where CH4 is prominent. Unfortunately, this pre- 
vents us from estimating the abundance of CH4 - 
and therefore of carbon - in Gl 229B. 

Ammonia shows significant absorption in both 
the H and K band spectra. The line list for NH3 
starts at A > 1.415 /tm. Because the line list does 
not include transitions from highly excited levels 
which occur at T > 300 K, the NH3 opacity we 
compute at a given wavelength is strictly a lower 
limit to the actual opacity. 

Except for the collision-induced absorption by 
H2, the most important molecular absorber in Gl 
229B is H 2 0, for which the opacity is now rela- 
tively well understood. We use the most recent 
and most extensive ab initio line list (Partridge & 
Schwenke (1997), 3 x 10 8 transitions). This line 
list is essentially complete for T ^ 3000 K. As a 
demonstration of the equality of this database, we 
find that the H2O features computed with this line 
list correspond extremely well in frequency with 
the observed features of Gl 229B (e.g. Figs. 7 to 
9). However, we find noticeable discrepancies in 
the relative strengths of H2O features which we 
attribute to the calculated oscillator strength of 
the transitions (§5.1). This effect can also be seen 
in Fig. lc of Griffith et al. (1998). At high reso- 
lution, the distribution of molecular transitions in 
frequency and strength is nearly random, and the 
inaccuracies in oscillator strengths we have found 
limit the accuracy of model fitting in a fashion 
similar to noise. This "opacity noise" is at least 
as significant as the noise intrinsic to our data. 

5. Results from fitting the spectra 

We have constructed a grid of synthetic spectra 
for the three models shown in Fig. 4 with metal- 
licity -0.7 < [M/H] < 0.1 in steps of 0.1. We use 
these modeled spectra to fit four distinct spectral 
regions (the "red", J, H, and K spectra) sepa- 
rately to determine the metallicity as a function 
of gravity with an internal precision of ±0.1 dex. 



5 



For the purpose of fitting the data, the synthetic 
spectra were renormalized to the observed flux at 
a selected wavelength in each spectral region. 

The model spectra show some distortions in the 
overall shape of the spectrum which are probably 
due to remaining uncertainties in the (P, T) profile 
of the atmosphere, the inadequate CH4 opacities, 
and possible effects of dust opacity. Considering 
the additional problems with the strength of the 
H 2 features, we elected to do all fits "by eye," 
except where otherwise noted. We discuss the fit- 
ting of each spectral interval below. In the interest 
of brevity, we present a detailed discussion of fits 
obtained only with the model of intermediate grav- 
ity (model B). The best fits obtained with models 
A and C are very nearly identical to those with 
model B. The results are summarized in Table 2. 

5.1. The "red" spectrum 

The "red" spectrum extends from 0.83 to ~ 
f and is formed deep in the atmosphere where 
1100 < T < 1500 K. The spectra of Schultz ct al. 
(1998) and Oppcnhcimer ct al. (1998) reveal two 
lines of Cs I (at 0.852 and 0.894 /zm) and a strong 
H 2 band but not the bands of TiO and VO com- 
mon to late M dwarfs and early L dwarfs (Fig. 
7). Refractory elements, such as Ti, Fe, V, Ca, 
and Cr, are expected to be bound in condensed 
compounds in a low-temperature atmosphere such 
as that of Gl 229B and therefore are not avail- 
able to form molecular bands (Feglcy & Lodders 
1996; Marley et al. 1996; Burrows & Sharp 1999). 
With the exception of the strong, unidentified fea- 
ture at 0.9874 ^m, all features between 0.89 and 
1.0 nm can be attributed to H2O. An overlap of 
a weak band of H 2 and a weaker CH4 band 
causes the small depression at 0.894 fim. which was 
tentatively tentatively attributed to CH4 by Op- 
penheimer et al. (1998) and Schultz et al. (1998). 
Features below 0.89 /im cannot be identified at 
present. The two Cs I lines are not included in 
our model. 

The flux from Gl 229B is also observed to de- 
crease very rapidly toward shorter wavelengths 
(Schultz et al. 1998; Oppenheimer et al. 1998; 
Golimowski et al. 1998), which, in the absence of 
the strong TiO and VO bands, is evidently caused 
by the presence of a missing source of opacity in 
the atmosphere. Spectra computed with molecu- 
lar opacities only (but excluding TiO, VO, FcH, 



etc) predict visible fluxes which are grossly over- 
estimated (Griffith ct al. 1998) but the detailed se- 
quence of absorption features of the spectrum are 
well reproduced, indicating that the short wave- 
length flux is suppressed by a smooth opacity 
source. 

We fit the red spectrum of Gl 229B between 
0.82 and 1.15 /xm to obtain the metallicity. We 
have recalibrated the published spectrum (Oppen- 
heimer et al. 1998; Geballe et al. 1996; Leggett et 
al. 1999) using the HST photometry of Golimowski 
et al. (1998). We model the continuum opac- 
ity with a layer of condensates following the ap- 
proach of Griffith et al. (1998). Condensates are 
expected in the atmosphere of Gl 229B on the ba- 
sis of chemical equilibrium calculations (Lodders 
1999; Burrows & Sharp 1999) and can provide the 
required opacity. Alternatively, Tsuji, Ohnaka & 
Aoki (1999) and Burrows, Marley, & Sharp (1999) 
attribute this rapid decline to the pressure broad- 
ened red wing of the 0.77 /xm K I resonance dou- 
blet. The first optical spectrum of a T dwarf 
(SDSS 1624+0029) shows that the latter expla- 
nation is correct (Liebert et al. 2000). The nature 
of this opacity source is not very important for 
the determination of the metallicity, however, as 
long as the proper continuum opacity background 
is present in the calculation. 

The dust opacity is computed with the Mie the- 
ory of scattering by spherical particles and is de- 
termined by the vertical distribution of the parti- 
cles, their grain size distribution, and the complex 
index of refraction of the condensate. The cloud 
model of Griffith et al. (1998) is described by 1) 
the vertical density profile of the condensate, taken 
as: 

n d = AP, 

where n d is the number density of condensed par- 
ticles, P is the ambient gas pressure, and the cloud 
layer is bound by 

■Ptop < P < 100 bar; 

2) the size distribution of the particles 



f(d) = ^ exp 



Hd/dp) 
lntr 



where d is the diameter of the particles; and 3) the 
complex index of refraction of the condensate 

n(A) = 1.65 + m;(A). 
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The parameters Ptop, A, do, o~, the function rii(A), 
and the metallicity of the atmosphere are free pa- 
rameters. 

Such a multi-parameter fit of the observed spec- 
trum is not unique. Furthermore, arbitrarily good 
fits of the "continuum" flux level can be obtained 
by adjusting the imaginary part of the index of re- 
fraction since its wavelength dependence is weakly 
constrained a priori. Our results for the three 
models are qualitatively similar to those of Griffith 
ct al. (1998). Typical values of the fitted dust pa- 
rameters are A = 415cm~ 3 bar~ 1 , P top = 0.5 bar, 
do = 0.21 /xm, and a — 1.3, with the imaginary 
part of the index of refraction decreasing from 
n t = 1.70 at 0.8 /xm to 0.01 at 1.10 /xm. These 
parameters applied to model B with a metallicity 
of [M/H]=— 0.3 result in the fit shown in Figure 7. 

The metallicity of the atmosphere [M/H] is 
largely independent of the dust parameters, how- 
ever, as it is constrained by the amplitude of the 
features in the H 2 band which we fit between 
0.925 and 0.98 /an. A larger metallicity results in 
a larger amplitude of the features inside the band. 
We obtain the same value of [M/H] as long as a 
good fit of the "continuum" flux level is obtained, 
regardless of the particular values of dust param- 
eters. The lower panel of Fig. 7 clearly shows 
differences in the relative strengths of the absorp- 
tion features in the H2O band between the syn- 
thetic and observed spectra. Similar differences 
also occur in the J, H, and K bands. The same 
differences are found for all three atmospheric pro- 
files and point to inaccuracies in the oscillator 
strength of the ab initio line list of H2O (Partridge 
& Schwenke 1997). The metallicity is fitted to 
give the best overall fit of these features with a 
precision of ±0.1 dcx. 

Strictly speaking, this procedure gives the H2O 
abundance, or [O/H] rather than [M/H]. For solar 
metallicity, the condensation of silicates deep in 
the atmosphere of Gl 229B will reduce the amount 
of oxygen available to form H 2 by ~ 15% (Fegley 
& Prinn 1988) which implies that 

[M/H] = [O/H] + 0.07. (1) 

The correction, which decreases for subsolar 
metallicities, is smaller than our fitting uncer- 
tainty and will be ignored hereafter. 



5.2. J band spectrum 

The J band spectrum probes the most trans- 
parent window of the spectrum of Gl 229B and is 
formed at great depth where 1500 < T < 1600 K 
and P ~ 30 bar (Fig. 6). Our spectrum contains 
almost exclusively H 2 features, with the excep- 
tion of CH4 absorption for A Ss 1.215 /xm and of 
two prominent K I lines (Fig. 1). Since we have 
elected to ignore CH4 bands and our synthetic 
spectra do not include alkali metal lines, we fit the 
J spectrum between 1.215 and 1.298 /xm to deter- 
mine the metallicity from the depth of the H 2 
absorption features. Figure 8 shows the effect of 
the metallicity on the spectrum (top panel) and 
our best fit (bottom panel) for model B. The flux 
level in the J spectrum varies by a factor of 3 and 
our best fit shows distortions in the general shape 
of the spectrum. The distortions may be caused by 
a combination of uncertainties in the (T, P) profile 
of the atmosphere, problems with the H 2 opac- 
ities or a small amount of dust opacity in the in- 
frared spectrum of Gl 229B (not modeled). Since 
the fit is based on the depth of the features, we ig- 
nore these distortions and fit the logarithm of the 
flux rather than the flux, as shown in Fig. 8. As in 
the red spectrum (§5.1), we find a remarkable cor- 
respondence of spectral features between the ob- 
served and modeled spectra although the model 
spectrum is somewhat less successful at reproduc- 
ing the relative strengths of the H 2 features. 

5.3. H band spectrum 

While the H band spectrum falls between the 
red side of a strong H 2 band (for A & 1.59 /xm) 
and the blue side of a prominent CH 4 band (for 
A 1.59 /im), NH3 and H 2 S have non-negligible 
opacity in this wavelength interval which also in- 
cludes a band of CO. 

The strong CH4 band is responsible for the 
turnover of the flux at 1.59 /tm. In this wavelength 
range, the CH4 opacity is described in our calcu- 
lation by the Strong et al. (1993) laboratory mea- 
surements which are restricted to T = 300 K. As 
a consequence, the CH4 band comes in at 1.61 /tm 
in the synthetic spectra, which results in strong 
departures of the mean flux level between data 
and models for A > 1.58 /xm. This limits our 
analysis of the H band spectrum to wavelengths 
shorter than 1.58 /xm. In this wavelength range, 
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the spectrum is formed deep in the atmosphere, 
where 1200 < T < 1350 K and P ~ 10 bar (Fig. 
6). Within this spectral region, we determine the 
abundance of NH3 as a function of the mctallicity, 
but cannot untangle the two. We also comment 
on the presence of H 2 S and CO. 

5.3.1. Metallicity and ammonia 

While NH 3 absorption can significantly affect 
the slope of the spectrum shortward of 1.56 /xm 
(Fig. 9) , there is no distinctive feature at this spec- 
tral resolution (R = 2100) to provide an unam- 
biguous detection of this molecule. The detection 
of CO at 4.7 /an well above the equilibrium abun- 
dance (Noll et al. 1997) suggests the possibility 
that NH 3 may also depart from its chemical equi- 
librium abundance (Fegley & Loddcrs 1994). We 
therefore vary the abundance of NH3 by reducing 
its chemical equilibrium abundance - as computed 
for a given metallicity - uniformly throughout the 
atmosphere by constant factor. We found that re- 
duction factors of 1, 0.5, 0.25, and provide an 
adequate grid of NH3 abundances given the S/N 
ratio of the data and the residual problems with 
the H 2 opacity. 

The effect of varying the metallicity on the H 
band spectrum is shown in Fig. 10 for model B. All 
spectra in Figure 10 were computed with the equi- 
librium abundance of NH3 . In this case, we find a 
best fitting mctallicity of [M/H]= -0.3. Figures 9 
and 10 show that varying the abundance of NH 3 
and varying the metallicity have very similar ef- 
fects on the synthetic spectrum. In the absence 
of any distinctive feature of NH 3 , it is not possi- 
ble to determine both the metallicity and the NH3 
abundance separately. For each value of [M/H], we 
can adjust the NH 3 abundance to obtain a good 
fit, higher metallicities requiring lower NH3 abun- 
dances. The best fitting solutions are given in Ta- 
ble 2. These fits are nearly indistinguishable from 
each other although the higher metallicity fits are 
marginally better. 

It is in the H band that we find the poorest 
match in the detailed features of the data and the 
models. While the fit is determined by matching 
the relative amplitudes of the features, the two- 
parameter fit (mctallicity and NH 3 abundance) we 
have performed amounts to little more than fitting 
the slope of the spectrum. As an internal check on 
our fitting procedure and precision, we have have 



verified that our fits indeed have the same slope as 
the data by plotting data and models at a very low 
spectral resolution which eliminates all absorption 
features. 

5.3.2. Hydrogen sulfide 

Hydrogen sulfide (H 2 S) has non- negligible 
opacity over most of the wavelength range of our 
fit to the H band. Figure 11 shows two spec- 
tra computed with and without H 2 S opacity. Its 
opacity is weaker than that of NH3 however, and 
there are no distinctive features which would allow 
a positive identification. Since the H 2 S features 
are fairly uniformly distributed in strength and 
wavelength, fitting the spectrum with either the 
chemical equilibrium abundance of H 2 S or with 
no H 2 S at all only has a small-to- negligible effect 
on the determination of the NH3 abundance for 
a given metallicity (Table 2). For the low gravity 
model A, the NH3 abundance determined from 
spectra without H 2 S is about half of the value 
found with the chemical equilibrium abundance of 
H 2 S. The difference decreases at higher gravities 
and is negligible for model C. 

Chemical equilibrium calculations (Fegley & 
Lodders 1996; Lodders 1999) indicate that H 2 S 
is present in Gl 229B with an abundance essen- 
tially equal to the elemental abundance of sul- 
fur in the atmosphere (see §6.3 for further discus- 
sion). Unfortunately, we arc unable to ascertain 
the presence of H 2 S in the atmosphere of Gl 229B 
at present. 

5.3.3. Carbon monoxide 

The discovery of CO in the 4.7 yum spectrum 
of Gl 229B with an abundance about 3 orders of 
magnitude higher than predicted by chemical equi- 
librium revealed the importance of dynamical pro- 
cesses in its atmosphere (Fegley & Lodders 1996; 
Noll et al. 1997; Griffith and Yelle 1999). The 
4.7 /um spectrum probes the atmosphere at the 2 
- 3 bar level where T ~ 900 K (Fig. 6). At this 
level, chemical equilibrium calculations predict a 
CO abundance of A co = 4.7 x 10~ 8 while Noll 
et al. (1997) found 5 x 10~ 5 <, X co ^ 2 x 10~ 4 . 
The second overtone band of CO falls within the 
H band and, in principle, could provide a deter- 
mination of Xqo at a deeper level of the atmo- 
sphere, where P — 14 bar and T ~ 1400 K (Fig. 
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6). Figure 12 shows a comparison of the data with 
synthetic spectra computed with various amounts 
of CO for model B with solar metallicity. The first 
two spectra are computed with the chemical equi- 
librium abundance of CO (X C o = 4.90 x at 
14 bar) and without CO (Xco = 0). These two 
spectra arc very nearly identical. A third spec- 
trum is computed in the unrealistic limit where all 
the carbon in the atmosphere is in the form of CO 
(Xco = 2.97 x 1CP 4 ), which represents the maxi- 
mum possible CO enhancement. As we found for 
NH 3 and H 2 S, there is no distinctive spectral sig- 
nature of CO at this resolution (R = 2100). Our 
extreme case represents a flux reduction of ~ 2a 
at best. We are unable to constrain the CO abun- 
dance with our data. Since the iS.v = 3 band of 
CO is <~ 10 3 times weaker than the fundamental 
band at 4.7 /im, obtaining a useful CO abundance 
from H band spectroscopy will be a difficult un- 
dertaking. 

5.4. K band spectrum 

Of our three near infrared spectra, the K band 
spectrum is formed highest in the atmosphere: 
T = 950 - 1020 K and P ~ 3 bar (Fig. 6). This 
is the same level as is probed with 4.7 /im spec- 
troscopy. As in the H band, the spectrum contains 
mainly H 2 features on the blue side (A ^ 2.1 /im) 
and a strong CH4 band appears at longer wave- 
lengths. There are also several NH 3 features for 
A < 2.05 ^im and the models predict an isolated 
feature of H 2 S (Fig. 3). 

Figure 13 compares a spectrum computed for 
model B with [M/H]= -0.3 with the entire K 
band spectrum. There is an excellent agreement 
in the structure of the spectrum even though the 
overall shape is not very well reproduced. For 
A < 2.05 /im, the model predicts strong features 
of NH 3 which we discuss in the next section. Be- 
yond 2.12 ^m is a CH4 band which is too weak 
in the model. The structure within the modeled 
band is remarkably similar to the observed spec- 
trum, however. This much better agreement of the 
CH4 band than we obtained in the H band is due 
to two factors. First, in this band the CH4 opac- 
ity is computed from a line list, and can therefore 
be computed as a function of temperature rather 
than at a fixed value of 300 K. Second, the lower 
temperature where the band is formed (Fig. 6) re- 
duces the effect of the incompleteness of the CH 4 



line list above 300 K. Nevertheless, we do not con- 
sider the CH 4 features here and limit our analysis 
to A < 2.10 /im. 

The K band spectrum provides a unique op- 
portunity: Once the metallicity is determined by 
matching the depth of the H2O features between 
2.05 and 2.10 /im, the abundance of NH3 can be 
obtained by fitting its features below 2.05 /im. The 
fit of the metallicity is shown in Fig. 14 for model 
B, which shows [M/H]=0 and —0.5 (top panel) and 
our best fit, [M/H] =— 0.3 (bottom panel). Values 
for models A and C are given in Table 2. All fea- 
tures in this region are due to H2O and, as we 
found in the red spectrum and in the J and H 
bands, the oscillator strengths of the H 2 line list 
do not reproduce the relative strength of the ob- 
served features very well. 

5.4-1- Ammonia 

In the wavelength range shown in Fig. 15, the 
spectrum consists of a few NH3 features on a back- 
ground of H 2 absorption. Synthetic spectra pre- 
dict seven strong NH3 features in this spectrum, 
three of which are clearly present (2.033, 2.037 and 
2.046 /im), one is absent (2.041 /im) and three ap- 
pear to be missing (2.026, 2.029 and 2.031 /im). 
This constitutes an ambiguous detection of NH 3 
in Gl 229B. The determination of the abundance 
of NH3 from the K band spectrum is hampered by 
the limited accuracy of the oscillator strengths of 
the H 2 line list and by the incompleteness of our 
NH3 line list for temperatures above 300 K. The ef- 
fect of the former can be seen in the trio of features 
at 2.026, 2.029 and 2.031 /im, which overlap H 2 
absorption features. Even after removing all NH3 
opacity, these features are still too strong in the 
calculated spectrum (top panel of Fig. 15). The 
top panel of Fig. 15 as well as Fig. 13 show that 
for model B, the abundance of NH 3 derived from 
the chemical equilibrium for the adopted metal- 
licity of [M/H]= —0.3 is too high. We therefore 
consider a depletion in NH 3 in the atmosphere 
of Gl 229B at the level probed by the K band 
spectrum. Following the approach used in fitting 
the H band spectrum, we express this depletion 
as a fraction of the chemical equilibrium abun- 
dance of NH3 for the metallicity obtained inde- 
pendently from the amplitude of the H 2 features 
in the K band. This fraction is applied uniformly 
throughout the atmosphere for the computation 
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of the synthetic spectrum. Given the ambiguous 
presence of NH3 in the K band, we have deter- 
mined the optimal NH3 abundance by minimizing 
the x 2 of the spectral fit for 2.022 < A < 2.049 fxm. 
This gives a depletion factor of ^ 0.4 with no NH 3 
present being an acceptable fit. Restricting the 
fit to the region where NH3 features are clearly 
observed (2.032 < A < 2.049) gives a similar re- 
sult but favors a finite value for the NH 3 deple- 
tion of ~ 0.2. The results are summarized in Ta- 
ble 2. The lower panel of Figure 15 shows the 
model B fit obtained by reducing the NH 3 abun- 
dance throughout the atmosphere to ~ 25% of its 
chemical equilibrium value for the adopted metal- 
licity. With this significant degree of depletion, 
the model reproduces the three detected features 
(2.033, 2.037 and 2.046 ^m) extremely well and 
makes the 2.041 /im feature consistent with the ob- 
servations. 

Because our NH 3 line list is incomplete at high 
temperatures, the opacity which we compute is 
strictly a lower limit to the actual NH 3 opacity at 
any wavelength. If the incompleteness is signifi- 
cant for the features found in the K band spec- 
trum, then the NH 3 abundance is actually lower 
still. It appears extremely unlikely that the er- 
rors in the oscillator strength of the H 2 tran- 
sitions would conspire to mimic the depletion of 
NH 3 which we find. For example, if we imagine 
that there is no NH 3 depletion (dotted curve in 
the top panel of Fig. 14), the residuals between 
the data and the fitted spectrum for A < 2.05 /im 
would be much larger than the typical mismatch 
which we find in H 2 features in all four spectra 
presented here. We consider the depletion of NH 3 
in the K band spectrum to be firmly established. 

5.4-2. Hydrogen sulfide 

We are not able to establish the presence of H 2 S 
from our H band spectrum (§5.3.2). Throughout 
the K band, the H 2 S opacity is generally over- 
whelmed by H 2 and CH4 absorption. However, 
there is a peak in the opacity of H 2 S which is about 
one order of magnitude higher than all other opac- 
ity maxima (Fig. 16). Our synthetic spectra indi- 
cate that this feature is strong enough to become 
visible in the midst of the background of H 2 and 
CH4 features. Figure 17 shows the relevant por- 
tion of the K band spectrum with synthetic spec- 
tra for all three models (A, B, and C) using the 



fitted metallicity (Table 2). All three panels are 
remarkably similar. The strength of the predicted 
feature is well above the noise level of the data, and 
taken at face value, Figure 17 indicates a proba- 
ble depletion of H 2 S by more than a factor of 2. 
On the other hand, chemical equilibrium calcu- 
lations indicate that the H 2 S abundance should 
be very near the elemental abundance of sulfur, 
even in the presence of vertical transport and con- 
densation (Fegley & Lodders 1996). Since there 
is no reason to expect a significant depletion of 
H 2 S, the discrepancy is probably due to remain- 
ing uncertainties in the opacities. Our H 2 S line 
list is based on an ab initio calculation (R. Watt- 
son, priv. comm.) which hasn't been compared to 
laboratory data in this part of the spectrum. The 
strong feature centered at 2.1084 /im is a blend 
of three strong lines from three different bands of 
H 2 S. Possible errors in the position or strength of 
these lines could significantly reduce the ampli- 
tude of the feature in our synthetic spectra. The 
limitations of the background opacity of H 2 and 
CH4 may also be responsible for the observed mis- 
match. Nevertheless, it is desirable to look for this 
feature at a higher resolution and a higher S/N ra- 
tio as an absence of sulfur in Gl 229B would be a 
most intriguing result. 

5.5. Carbon monoxide at 4.7/im 

Given our determination of T c g and [M/H] in 
terms of the surface gravity, we can obtain the 
abundance of CO from the 4-5 /xm spectra of Noll 
ct al. (1997) and Oppenheimer et al. (1998) which 
is consistent with our results. For each model - A, 
B, and C - and using the metallicity given in Ta- 
ble 2, we computed synthetic spectra with various 
CO abundances. The latter is varied freely with- 
out imposing stoichiometric constraints. The syn- 
thetic spectra are binned to the wavelength grid 
of the data and fitted to the data by a normal- 
ization factor adjusted to minimize the x 2 - The 
X 2 of the fitted spectra shows a well-defined min- 
imum as a function of the CO abundance, Aco 
(Table 3). The uncertainty on the CO abundance 
is obtained by generating synthetic data sets by 
adding a gaussian distribution of the observed 
noise to the best fitting model spectrum. After 
doing the same analysis on 1000 synthetic data 
sets, we obtain a (non-Gaussian) distribution of 
values of Aco- The uncertainties given in Table 3 
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correspond to the 68% confidence level. The Op- 
penheimer et al. (1998) spectrum gives CO abun- 
dances which are 0.1 dex higher than those ob- 
tained from the Noll et al. (1997) spectrum; which 
is well within the fitting uncertainty. Noll et al. 
(1997) found a CO mole fraction of 50 to 200 ppm 
(—4.3 < logXco < —3.7) by assuming a H 2 
abundance of 300 ppm (effectively, [M/H]= —0.6), 
which agrees well with our result for model A 
which has [M/H]= —0.5. Our results are also 
consistent with those of Griffith and Yelle (1999) 
who find CO abundances of > 20 and > 100 ppm 
(logX C o > -4.7 and -4) for [M/H]=-0.6 and 0, 
respectively, using the data of Noll et al. (1997). 

Since this part of the spectrum contains only 
H 2 and CO features, our fitting procedure is sen- 
sitive only to the CO to H 2 abundance ratio. 
Chemical equilibrium calculations show that the 
H2O abundance scales linearly with the metallicity 
at the level probed with 4-5 /mi spectroscopy (i.e. 
all the gas phase oxygen is in H2O). Accordingly, 
the CO abundance we find scales with the metal- 
licity of the model. As shown in Figure 18, the 
CO abundance determined from the 4.7 /mi band 
corresponds approximately to the CO/CH4 tran- 
sition in chemical equilibrium. The observations 
definitely exclude the very-low chemical equilib- 
rium abundance of CO. The fact that the extreme 
case where all carbon is CO provides an acceptable 
fit to the data (while we know that a good fraction 
of the carbon is in CH 4 ) is due to the rather noisy 
spectra. 

6. Discussion 

6.1. Surface Gravity 

It is highly desirable to restrict the surface grav- 
ity g of Gl 229B to an astrophysically useful range. 
Since Lboi is known, a determination of g fixes T e g, 
the radius, the mass, and the age of Gl 229B (Fig. 
4), as well as the metallicity and the abundances 
of important molecules such as CO and NH 3 . The 
large uncertainty on g results in a large uncer- 
tainty in the mass of Gl 229B and on the age of the 
system determined from cooling tracks. While a 
dynamical determination of the mass may be pos- 
sible in a decade or so (Golimowski et al. 1998), 
a spectroscopic determination might be obtained 
much sooner. 

Unfortunately, it is not possible to constrain the 



gravity better than logg = 5 ± 0.5 with the data 
and models presently available. Our high resolu- 
tion spectroscopy does not allow us to choose be- 
tween models A, B and C (Table 2) as an increase 
in gravity can be compensated by an increase in 
metallicity to lead to an identical fit. The grav- 
ity sensitivity of the K band synthetic spectrum 
models reported by Burrows et al. (1999a) occurs 
for a fixed metallicity only. 

Similarly, Allard et al. (1996) and Burrows et 
al. (1999a) report that the spectral energy distri- 
bution of Gl 229B models is fairly sensitive to the 
gravity. But this is true only for a fixed metallic- 
ity. For the three models indicated in Fig. 4, and 
using the metallicity we have determined for each 
(Table 2) the gravity dependence of the infrared 
colors is A(H - K)/A(\ogg) = 0.03 and A(J - 
K)/A(logg) = 0.03 and A(K - L')/A(logg) = 
—0.08 (Table 4). This dependence is very weak 
in the light of the uncertainty in the photometry 
of Gl 229B. More problematic is the fact that the 
synthetic H — K disagrees with the photometry. 
Furthermore, the incomplete CH4 opacities used 
in the spectrum calculation almost certainly result 
in an inaccurate redistribution of the flux in the 
near infrared opacity windows which determines 
the broad band colors. An example of this effect 
on the H band flux can be seen in Fig. 10. We 
conclude that a photometric determination of the 
gravity is not possible at present. 

An alternative approach to determining the 
gravity of Gl 229B is through a study of the 
pressure-broadened molecular lines of its spec- 
trum. The spectrum of Gl 229B is formed of a 
forest of unresolved molecular lines maily due to 
H2O, CH4, and NH3. Because of the limitations of 
the CH4 and NH3 opacity data bases, a detailed 
study of molecular features is best performed in 
spectral domains where these two molecules do not 
contribute significantly to the opacity. Spectro- 
scopic observations with a resolution of ~ 50 000 
can reveal the shape of individual H2O lines in re- 
gions where they are relatively sparse, e.g. from 
2.08 to 2.105 /mi. 

6.2. Metallicity 

Our determination of the metallicity of Gl 
229B, with an uncertainty of ±0.1, is given in 
Table 2 for the three gravities considered. We find 
an excellent agreement between our three indc- 
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pendent determinations of [M/H] for each gravity 
and conclude that Gl 229B is likely to be depleted 
in heavy elements, e.g. oxygen. The metallicity 
is near solar at high gravity and decreases signifi- 
cantly for lower gravities. In their analysis of the 
"red" spectrum, Griffith et al. (1998) found a H 2 
abundance between 0.3 and 0.45 of the solar value 
for a T e ff = 900 K, log g = 5 model and they adopt 
a value of 0.25 in Griffith and Ycllc (1999). 

Metal depletion in Gl 229B is consistent with 
the analysis of the 0.985 - 1.02 /xm FeH band in the 
spectrum of the primary star Gl 229A by Schiavon, 
Barbuy & Singh (1997) who find [Fc/H]= -0.2. 
The infrared colors of Gl 229A imply that it is 
slightly metal- rich, however (Leggett 1992). The 
relative metallicity of the components of this bi- 
nary system may have been affected by their for- 
mation process. If the pair formed from the frag- 
mentation of a collapsing cloud (like a binary star 
system), the two objects should share the same 
composition. If the brown dwarf formed like a 
planet, from accretion within a dissipative Keple- 
rian disk around the primary, the selective accre- 
tion of solid phase material could lead to an en- 
richment in heavy elements compared to the pri- 
mary star, as is observed in the gaseous planets 
of the solar system. The low mass of the primary 
(~ 0.5M Q ) and the large semi-major axis and ec- 
centricity (a ^ 32 AU and e ^ 0.25, Golimowski 
et al. (1998)) suggests that the binary formation 
process, and therefore equal metallicities, are more 
plausible. A more detailed study of the metallicity 
of Gl 229A is desirable to better understand the 
history of this system. 

6.3. Atmospheric chemistry and molecu- 
lar abundances 

The results of our analysis of the metallicity 
and abundances of several molecules in the atmo- 
sphere of Gl 229B are summarized in Fig. 18. 
Each panel corresponds to a different model (see 
Table 2) and shows the abundances of important 
molecules as a function of depth in the atmosphere 
based on chemical equilibrium calculations includ- 
ing condensation cloud formation. The chemistry 
of these abundant molecules is fairly simple. The 
abundance of H2O is uniformly reduced by ~ 16% 
by silicate condensation. Except for a small deple- 
tion for log T ^ 3 due to the condensation of Na 2 S, 
all sulfur is found in H 2 S. The other molecules 



shown are not affected by condensation. Nitrogen 
is partitioned between N 2 and NH 3 , with the lat- 
ter being favored at lower temperatures and higher 
pressures. NH3 dominates near the surface and 
rapidly transforms into N 2 at the higher temper- 
atures found deeper in the atmosphere. Deep in 
the atmosphere, the higher pressures cause a par- 
tial recombination of NH3 and the ratio of the NH3 
to N 2 abundances increases slowly with depth. In 
a similar fashion, all elemental carbon is found in 
CH4 at the surface but CO starts to form at higher 
temperatures and rapidly becomes the most abun- 
dant carbon-bearing molecule. The formation of 
CO consumes H 2 0, as can be seen in Fig. 18. 

In each panel, a dotted box indicates the CO 
abundance we have determined from the 4.7 /xm 
spectrum, the location of the box along the or- 
dinate shows the level probed at this wavelength 
(Fig. 6). As originally discussed by Noll et al. 
(1997) and Griffith and Ycllc (1999), the CO abun- 
dance is ~ 3 orders of magnitude larger than the 
chemical equilibrium value. Stochiomctric con- 
straints imply that this also results in a signifi- 
cant reduction of the CH4 abundance at the 870 
- 950 K level in Gl 229B. In model B, the equilib- 
rium abundance of CH 4 at 900 K is 2.97 x 10~ 4 , 
which corresponds to the abundance of elemen- 
tal carbon. The CO abundance determined from 
the 4.7 //m band is 1.6 x 10~ 4 (with a large error 
bar). Conservation of the total number of carbon 
atoms then requires that the CH4 abundance be 
<~ 1.4 x 10~ 4 , a full factor of 2 below the equilib- 
rium abundance. If we use the lowest CO abun- 
dance allowed by our analysis, a 25% reduction 
of CH 4 relative to its equilibrium abundance at 
900 K results. Depletion of CH 4 at this depth is 
readily accessible spectroscopically in the 1.6 and 
3.3 ^m bands, and may also affect the 2.3 ^im band 
if the non-equilibrium CO abundance persists at 
higher levels (Fig. 6). After H 2 0, CH 4 is the 
most important near infrared molecular absorber 
in Gl 229B. Accurate modeling of the spectrum de- 
mands a careful treatment of the non-equilibrium 
CH4 abundance. 

Similarly, solid boxes show the NH3 abundance 
determined from our H and K band spectra. 
While the H band abundance is in excellent agree- 
ment with the equilibrium value, there is a clear 
depletion of NH3 in the K band. The H and K 
band abundances are marginally consistent with 
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each other but it appears that the NH3 abundance 
decreases upwards through the atmosphere. 

6.4. Non-equilibrium processes 

Processes which take place faster than the time 
scale of key chemical reactions can drive the com- 
position of the mixture away from equilibrium. 
The case of CO/CH4 chemistry has been well- 
studied in the atmosphere of Jupiter where an 
overabundance of CO is also observed. Carbon 
monoxyde is a strongly bound molecule and the 
conversion of CO to CH4 through the (schematic) 
reaction 

CO + 3H 2 -» CH 4 + H 2 (2) 

proceeds relatively slowly, while the reverse reac- 
tion is much faster. Vertical transport, if vigorous 
enough, can carry CO-rich gas from deeper levels 
upwards faster than the CO to CH4 reaction can 
take place. The CO/CH4 ratio at any level is fixed 
( "quenched" ) by the condition 

7"chem — 7"mix; (3) 

where T c h C m is the chemical reaction time scale and 
T m i x the dynamical transport time scale. Wher- 
ever r m i x < Tchom in the presence of a verti- 
cal gradient in the equilibrium abundance, non- 
equilibrium abundances will result. 

As discussed by Fegley & Lodders (1996), Noll 
ct al. (1997), Griffith and Yclle (1999) (and ref- 
erences therein), this naturally explains the very 
high CO abundance observed at the 900 K level. In 
this picture, CO-rich gas would be carried upward 
from T Jb 1400 K. Convection is the most obvious 
form of vertical transport in a stellar atmosphere 
but in Gl 229B the convection zone remains about 
3 pressure scale heights below the level where CO 
is observed (shaded area in Fig. 18). Perhaps 
convective overshooting can transport CO to the 
observed level. Griffith and Yelle (1999) propose 
"eddy diffusion" as a slower, yet adequate trans- 
port mechanism. The eddy diffusion (or mixing) 
time scale is constrained by the poorly known CO 
abundance and the somewhat uncertain chemical 
pathway between CO and CH4. From Griffith and 
Yelle (1999), we infer that T m ; x ^ 1 to 10 years and 
could be much smaller. 

In analogy to the CO/CH4 equilibrium, a low 
NH3 abundance can be explained by vertical 



transport which can quench the NH3/N2 ratio at 
a value found in deeper layers in the atmosphere. 
As it is carried upward, N 2 is converted to NH 3 
by the reaction 

N 2 + 3H 2 -» 2NH 3 (4) 

The N 2 molecule is very strongly bound, however, 
and this reaction proceeds extremely slowly at low 
temperatures, much more slowly than reaction (2). 
Thcrmochemical kinetic calculations of the chem- 
ical lifetime for conversion of N 2 to NH3 were per- 
formed as described in Fegley & Lodders (1994). 
The time scale for reaction (4) along the (P,T) 
profiles of models A, B, and C and for P > 1 bar 
is given by 

3 75 x 10 4 

logr chem = — — + 0.22 log g - 23.08, (5) 

where r c h e m is in year, T is the temperature in K, 
and g is the surface gravity in cm/s 2 . This time 
scale assumes that the N 2 conversion occurs in 
the gas phase although it could possibly be short- 
ened by catalysis on the surface of grains. The 
time scale increases very steeply with decreasing 
temperature. At the level probed in the H band, 
logT = 3.1 and logr c h C m = 7.8 and in the K band, 
logT = 2.95 and logT c hom = 20.1! At some in- 
termediate level (which depends on g), the time 
scale for the conversion of N 2 into NH3 becomes 
longer than the age of Gl 229B. In view of the rel- 
atively very short mixing time scale inferred from 
the CO abundance, it follows that the NH 3 abun- 
dance in the H and K bands is entirely determined 
by non-equilibrium processes, not by reaction (4) . 
At depths where logT > 3.23 (which corresponds 
to the top of the convection zone), T c h C m < 1 yr 
and the reaction proceeds fast enough to estab- 
lish chemical equilibrium between NH3 and N 2 . 
We therefore expect that the N 2 /NH 3 ratio will 
be quenched at its value at logT w 3.23, through- 
out the remainder of the atmosphere. While the 
time scale for eddy diffusion increases in the up- 
per levels of the atmosphere, the extremely long 
Tchem ensures that the NH 3 /N 2 ratio remains un- 
changed, regardless of how slowly the vertical mix- 
ing proceeds. For model B, this corresponds to 
logYNH 3 = —5.2, in perfect agreement with the 
abundance found in the H band (Fig. 18b). 

While the abundance of NH 3 determined from 
the K band spectrum is not very precise, it is 
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marginally consistent with vertical mixing. Fig- 
ure 18 shows that it is more likely that the K 
band abundance is smaller than found in the H 
band, however. On the other hand, the abun- 
dances shown in Fig. 18 and Table 2 are depletion 
factors which were applied uniformly to the chem- 
ical equilibrium abundance profile of NH3, which 
has a large vertical gradient. For consistency with 
the mixed atmosphere picture, we have therefore 
redetermined NH3 abundances using a constant 
abundance throughout the atmosphere and found 
logX N H 3 < -5.0 and -4.7 ±0.15 from the K and 
H band spectra, respectively (model B). The for- 
mer is in good agreement with our simple predic- 
tion while the H band value is now rather high. 
The discrepancy between the H and K band re- 
sults thus persists in this new analysis. Changes 
in [M/H] within the ±0.1 uncertainty have little 
effect either. 

We consider the possibility that this vertical 
gradient in the NH 3 abundance may be caused 
by a different non-equilibrium process such as the 
photolysis of NH 3 by the UV flux from the primary 
star. Ammonia is a relatively fragile molecule 
which is easily dissociated by UV photons: 

NH 3 + hv -> NH 2 + H. (6) 

With a photodissociation cross section of « 6 x 
10~ 18 cm 2 per molecule, optical depth unity for 
the photodissociation of NH 3 is reached at pres- 
sures of a few millibars in Gl 229B. Photodisso- 
ciation of the much more abundant H2 molecules 
docs not effectively shield NH 3 from incoming UV 
photons since the two molecules absorb over differ- 
ent wavelength ranges. Photodissociation of NH 3 
therefore represents a net sink of NH 3 which oc- 
curs at the very top of the atmosphere. We can 
estimate a lower limit to the time scale of pho- 
todissociation of NH 3 by assuming that each pho- 
todissociating photon results in the destruction of 
a NH 3 molecule. The incident photon flux is 




(7) 

where photons causing dissociation of NH 3 are be- 
tween Ao and Ai, T\ is the flux at the surface of 
the primary star, R+ the radius of the primary 
star and d the separation of the binary system. For 
NH 3 , we have A ~ 1900 A and AA - 300 A (Moses, 



priv. comm.). The primary star has a dMl spec- 
tral type, with T\ ~ 10 7 erg cm _2 s _1 cm _1 and 
R* ~ 3.6 x 10 10 cm (Leggett et al. 1996). The bi- 
nary separation is d il 44 AU (Golimowski et al. 
1998). This results in N VY £ 8.3 x 10 3 cm- 2 s _1 . 
Photodissociation will affect significantly the NH 3 
abundance when 

(8) 

iVuv 

where a is the column density of NH 3 . Because 
the incident flux of UV photons is fairly low, this 
condition is satisfied only at pressures below a few 
microbars, i.e. very high in the atmosphere. In 
the region of interest, photodissociation destroys 
a very small fraction of NH 3 during one mixing 
time scale and therefore has little effect on the 
abundance of NH 3 . Photolysis of NH 3 cannot ex- 
plain the relatively low NH 3 abundance we find in 
the K band spectrum. 

We believe that the difference between the H 
and K band determinations of the NH 3 abundance 
arise from the limitations of the NH 3 opacity data 
used for the calculation of the synthetic spectra. 
As discussed above, the incompleteness of the NH 3 
line list for T > 300 K results in upper limits for 
the NH 3 abundances obtained by fitting the data. 
Since this effect increases with temperature, we 
expect that the abundance determined from the 
H band spectrum is overestimated relative to the 
one obtained from the K band spectrum, which 
is what we observe. Until NH 3 opacities become 
available for T ~ 800 - 1200 K, we will not be able 
to quantify this effect. 

Figures 6 and 18 show that ammonia offers a 
third window of opportunity for a determination 
of its abundance. The region between 8.3 and 
14.4 /im is rich in strong NH 3 features, the two 
strongest being at 10.35 and 10.75 /im (Fig. 19). 
This spectral region probes a higher level in the at- 
mosphere (T ~ 500 — 800 K) where the hot bands 
of NH 3 which arc missing from opacity data bases 
are less problematic than in the H and K bands. 
In this spectral region, the spectrum is very sen- 
sitive to the NH 3 /H 2 ratio, especially for NH 3 
abundances below 25% of the equilibrium value 
(Fig. 19). This would allow a good determination 
of the degree of NH 3 depletion in the upper lev- 
els of the atmosphere. We anticipate that 10 /im 
spectroscopy should reveal a NH 3 abundance of 
^ 10% of its equilibrium value. 
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7. Conclusion 

With the availability of extensive photometric, 
astrometric, and spectroscopic data, our picture of 
the atmosphere of Gl 229B is gradually becoming 
more exotic and more complex. The initial dis- 
covery of CH4 in its spectrum set it appart and 
has prompted the creation of a new spectral class, 
the T dwarfs. H 2 0, CO, Cs I, and K I have also 
been detected. There is good evidence that the 
rapid decrease of the flux at visible wavelengths 
is caused by unprecedently broad lines of atomic 
alkali metals (Liebert et al. 2000). The presence 
of condensates may also play a role in shaping the 
spectrum of Gl 229B. 

It is unfortunate that the surface gravity of Gl 
229B remains poorly constrained. We have not 
been able to further restrict the allowed range with 
our new J, H, and K spectroscopy. As a result, 
all our results are expressed as a function of grav- 
ity. This is the most significant obstacle to further 
progress in elucidating the astrophysics of this T 
dwarf. The surface gravity can probably be deter- 
mined from the study of the pressure broadened 
shape of molecular lines. 

We have found good evidence for the presence 
of NH3 in the spectrum of Gl 229B, which was 
expected from chemical equilibrium calculations. 
We have been able to determine its abundance at 
two different levels in the atmosphere, and we find 
a significant deviation from chemical equilibrium. 
A similar situation has been found with CO pre- 
viously (Noll et al. 1997) and this abundance pat- 
tern can be explained by vertical mixing in the 
atmosphere. The extent of the convection zone 
is not sufficient to account for the abundances we 
find and the mixing may be due to overshooting 
or to less efficient eddy diffusion. We find that 
NH 3 photolysis is not important in shaping the 
spectrum of Gl 229B. Because NH3 can be ob- 
served in three different bands corresponding to 
three distinct depths in the atmosphere, an accu- 
rate determination of its abundance in each band 
provides information on the time scale of mixing 
as a function of depth. This is an unusual and 
powerful diagnostic tool which can provide valu- 
able clues for modeling the vertical distributions 
of possible condensates. In principle, any absorber 
with a large abundance gradient through the visi- 
ble part of the atmosphere can be used to infer the 



details of the mixing process. Among detected and 
abundant molecules, only CO and NH 3 satisfy this 
criterion. Chemical equilibrium calculations with 
rainout of condensates (Lodders 1999; Burrows et 
al. 1999b) show that we can expect significant ver- 
tical gradients in the abundances of atomic K, Rb, 
Cs, and Na as they become bound in molecules 
(KC1, RbCl, CsCl and Na 2 S, respectively) in the 
cooler, upper reaches of the atmosphere. Cesium 
and potassium have been detected in the spectrum 
of Gl 229B, and resonance doublets of K I and 
Na I appear to shape the visible spectrum. How- 
ever, the chemical timescales for alkali metals are 
so short that they should always remain in thermo- 
dynamic equilibrium (Lodders 1999). Therefore, 
they cannot serve as probes of vertical mixing in 
Gl 229B. 

Further progress in understanding the atmo- 
sphere of Gl 229B requires better opacities for CH4 
and NH 3 , and, to a lesser extent, of H 2 0. A more 
accurate determination of the CO abundance from 
4-5 fxm spectroscopy is very desirable and will 
require higher signal-to-noise spectroscopy than is 
currently available. Similarly, 10 fim spectroscopy 
to determine the NH 3 abundance for P & lbar, 
while difficult, is important. 

The issue of vertical mixing and departures 
from chemical equilibrium gains importance when 
we consider that the observed departure of CO 
from chemical equilibrium implies a significantly 
reduced CH4 abundance, by conservation of the 
abundance of elemental carbon. Similarly, our re- 
sults imply that NH 3 absorption in the 10 /<m re- 
gion is reduced. Because CH4 is a significant ab- 
sorber in the near infrared, as is NH 3 in the 10 (ira 
range, departures from equilibrium must be taken 
into account when accurate modeling of the atmo- 
sphere and spectrum of Gl 229B is desired. This 
new level of complexity compounds the exoticism 
and the challenges posed by T dwarfs. 

The astrophysics of Gl 229B is far richer than 
has been originally anticipated. Gl 229B is cur- 
rently the only T dwarf known to be in a binary 
system. There is no evidence that the illumina- 
tion from the primary star has a significant ef- 
fect on the state of its atmosphere and Gl 229B 
is most likely typical of isolated T dwarfs. It re- 
mains the brightest and by far the best studied 
of the seven T dwarfs currently known, but the 
list should expand to several dozens during the 
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next 2-3 years (Burgasser et al. 1999; Strauss 
ct al. 1999). The existing body of work on Gl 
229B points to the most rewarding observations to 
conduct on T dwarfs. The possibility of studying 
trends in the physics of T dwarf atmospheres as a 
function of effective temperature is a fascinating 
prospect. 
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Fig. 1. — J band spectrum. The spectral res- 
olution is R = 2400 and the noise level at the 
maximum flux level is indicated by the error bar. 
Nearly all features in this spectrum are caused 
by H 2 absorption, except at short wavelengths 
where the CH 4 opacity becomes significant, and 
for two strong K I lines. The flux calibration is 
approximate. 

Fig. 2. — Same as Fig. 1 for the H band. The 
features seen in the spectrum are caused by H 2 
for A ^ 1.59 zxm. At longer wavelengths, the rapid 
decrease in flux corresponds to the blue side of a 
strong CH4 absorption band. Both NH 3 and H 2 S 
have significant opacity in this spectral range but 
show no distinctive feature at this resolution. The 
flux calibration is approximate. 

Fig. 3. — Same as Fig. 1 for the K band. The 
features seen in the spectrum are caused by H 2 
for A < 2.11 /xm. At longer wavelengths, the 
rapid decrease in flux corresponds to the blue 
side of a strong CH4 absorption band. Several 
strong NH3 features are seen shortward of 2.06 /xm. 
Models predict a single H 2 S absorption feature at 
2.1084 /xm. The flux calibration is approximate. 

Fig. 4. Evolution of solar metallicity brown 
dwarfs and giant planets in the effective temper- 
ature - gravity plane. The heavy solid lines are 
cooling tracks for objects with masses of 0.075, 
0.07, 0.06, 0.05, 0.04, 0.03, 0.02, and 0.01 M , 
from top to bottom, respectively (Burrows et al. 
1997). Evolution proceeds from right to left and 
isochrones (dotted lines) are labeled with the age 
in Gyr. The band crossing the center of the figure 
is the locus of all models with the luminosity of 
Gl 229 B (L bo i = 6.2 ± 0.55 x 1O- 6 L ). Filled 
symbols correspond to the models listed in Table 
2 and used for the present analysis. 

Fig. 5. — Temperature-pressure profiles of the 
three atmosphere models used in this analysis. La- 
bels correspond to the models shown in Fig. 5 and 
Table 2. 

Fig. 6. — The temperature and pressure at the 
photosphere (t\ = 2/3) as a function of wave- 
length for model B (T cff = 940 K, logg = 5). The 
temperature plotted (left axis) is equivalent to a 
brightness temperature. Horizontal bars indicate 
the wavelength range of NH 3 (solid lines) and CO 



(dotted lines) bands. The short vertical bar indi- 
cates a H 2 S absorption feature at 2.1084 /xm (see 
§5.4.2). 

Fig. 7. — A spectrum computed with model B, 
[M/H]= —0.3, and dust opacity (thin line) is com- 
pared to the observed spectrum of Oppcnheimcr 
et al. (1998) (heavy line). Top panel: Entire fitted 
region. Bottom panel: details of the H 2 absorp- 
tion band at a resolution R = 2250. 

Fig. 8. — Top panel: Effect of the metallicity on 
the J band spectrum. Data are compared to two 
model B spectra with [M/H]= -0.4 and 0. Bot- 
tom panel: The best fit for model B is obtained for 
[M/H]= — 0.2. See text. The synthetic spectra are 
normalized to the observations at 1.2472 /xm and 
are plotted at the same resolution (R = 2400). 

Fig. 9. — Effect of the NH 3 abundance on the 
H band spectrum. The spectrum is computed 
for model B with [M/H]= -0.3 with the chem- 
ical equilibrium abundance of NH3 (dotted line) 
and without NH 3 opacity (solid line). The dot- 
ted curve represents our best fit for this value of 
[M/H]. The synthetic spectra are normalized to 
the observations at 1.5744 /xm and plotted at the 
same resolution (R = 2100). 

Fig. 10. — Top panel: Effect of the metallicity 
on the H band spectrum. Data are compared to 
two model B spectra with [M/H]= —0.4 and 0. 
Bottom panel: The best fit for model B is obtained 
for [M/H]= —0.3. In all cases, the NH 3 abundance 
is determined from the chemical equilibrium. See 
text. See Fig. 9 caption for additional details. 

Fig. 11. — Effect of the H 2 S abundance on the 
H band spectrum. The spectra are computed for 
model B with [M/H]= —0.3 with the chemical 
equilibrium abundance of H 2 S and without H 2 S 
opacity. See Fig. 9 caption for additional details. 

Fig. 12. — Carbon monoxide absorption band in 
the H band. The spectra are computed for model 
B with [M/H]= with the chemical equilibrium 
abundance of CO (solid line) and with the CO 
abundance equal to the elemental carbon abun- 
dance (dotted line). A spectrum computed with 
the CO abundance set to zero overlaps the equilib- 
rium CO curve. See Fig. 9 caption for additional 
details. 
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Fig. 13. — Synthetic spectrum for model B with 
[M/H]= —0.3 compared to the observed K band 
spectrum. The synthetic spectrum is normalized 
to the data at 2.0513 /zm and plotted at the same 
resolution (R = 2800). 

Fig. 14. — Top panel: Effect of the metallicity 
on the K band spectrum. Data are compared to 
two model B spectra with [M/H]= —0.5 and 0. 
Bottom panel: The best fit for model B is obtained 
for [M/H]= —0.3. All features in this figure are 
caused by H 2 0. See Fig. 13 caption for additional 
details. 

Fig. 15. — Portion of the K band spectrum which 
shows NH 3 features. The data are shown by the 
heavy solid line, with a representative error bar. 
Ammonia features are indicated by tick marks. 
Top panel: Synthetic spectrum computed from 
model B with [M/H]= -0.3 (dotted line) and after 
removing all NH3 opacity (thin solid line). Bot- 
tom panel: Same as above but the equilibrium 
NH 3 abundance has been uniformly reduced by a 
factor of 0.25 throughout the atmosphere. See Fig. 
13 caption for additional details. 

Fig. 16. — Opacity of H 2 S in the 2 . 1 1 /im region for 
T = 800 K and P = 1 bar. The spectral resolution 
is R = 4.8 x 10 5 . Synthetic spectra predict that the 
strong feature at 2.1084 /^m should be observable. 

Fig. 17. — Hydrogen sulfide feature predicted at 
2.1084 /im. The observed spectrum is shown by 
a heavy solid line. The thin lines in each panel 
show spectra computed for the indicated model 
(A, B, and C). Each spectrum is computed with 
the best-fitting metallicity for that model (Table 
2). In each panel, the upper curve is calculated 
without H2S opacity, and the lower curve is com- 
puted with the chemical equilibrium abundance 
of H2S for the chosen metallicity. For model B 
(middle panel), the effect of reducing the equilib- 
rium H 2 S abundance by factors of 0.5 and 0.25, 
respectively, is also shown. All synthetic spectra 
are normalized to the data at 2.1062 /im and the 
resolution is R = 2800. 

Fig. 18. — Equilibrium chemistry in the atmo- 
sphere of Gl 229B. The curves show the abundance 
of the important molecular species as a function 
of depth for each of the three gravities, computed 
with the metallicity we have determined for each 



(Table 2). Depth increases toward the right along 
the ordinate axis and is indicated by both tem- 
perature and pressure scales. The extent of the 
convection zone is shown by the shaded region. 
Horizontal bars at the top show the regions probed 
by spectroscopic observations of various molecular 
bands. The top row of solid bars shows three re- 
gions where NH 3 is detectable: 10 fim. region and 
the K and H bands, from left to right, respec- 
tively. Dotted bars show the CO bands at 4.7 /im 
and in the H band, respectively. The vertical tick 
shows the depth where the 2.1084/im feature of 
H2S is formed. The dotted box shows the CO 
abundance determined from the at 4.7 /im band 
spectra of Noll et al. (1997) and Oppenheimer et 
al. (1998). Our determinations of the NH3 abun- 
dance in the K an H bands are shown by ver- 
tical boxes, the height of the box indicating the 
estimated uncertainty. For the K band measure- 
ments, an arrow indicates that the NH3 abundance 
could be close to zero. The labeling of the panels 
corresponds to the entries in Table 2. 

Fig. 19. — Ammonia features in the 10 /im re- 
gion. The spectra are computed with model B and 
[M/H]= —0.4. Each curve corresponds to a dif- 
ferent depletion of NH3 relative to its equilibrium 
abundance, the abundance of all other compounds 
being kept the same. The equilibrium abundance 
of NH 3 has been multiplied by 0, 0.25, 0.5, and 1, 
from top to bottom, respectively. All spectra are 
normalized at 10 /im and the resolution is R — 200. 
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Table 1 
Observing Log 



UT Date 


wavelength 
W 


Resolving 
Power 


Integ. Time 
(minutes) 


mean airmass 
Gl 229B 


mean airmass 
BS 1849 


19980125 


2.10-2.18 


3200 


80 


1.37 


1.39 


19980125 


2.02-2.10 


3100 


93 


1.40 


1.35 


19980126 


1.53-1.61 


2350 


53 


1.46 


1.39 


19980126 


1.25-1.30 


3000 


53 


1.33 


1.32 


19980126 


1.20-1.25 


2900 


53 


1.41 


1.40 


19980127 


2.02-2.10 


3100 


93 


1.36 


1.32 



Table 2 
Optimal parameters 



Model 


log g 

(cgs) 


(K) 




[M/H] a 




NH 3 /{NH 3 } cq 




0.92 ^m 


J H 


K 


H h 


H c 


K 


A 


4.5 


870 


-0.5 


-0.5 -0.3 


-0.5 





0.25 


& 0.40 










-0.4 




0.25 


0.5 












-0.5* 




0.5 


1 




B 


5.0 


940 


-0.3 


-0.2 -0.1 


-0.3 





<, 0.25 


£ 0.40 










-0.2 




0.25 - 0.5 


0.5 












-0.3* 




> 0.5 


1 




C 


5.5 


1030 


-0.2 


0.1 0.0 


-0.1 


0.5 


0.5 


£ 0.40 










-0.1* 




1 


1 





a An asterisk (*) indicates the adopted metallicity. 

b H 2 S abundance set to 0. 

C H 2 S abundance from chemical equilibrium. 



Table 3 

Abundance of CO from the 4.7 /zm band 



Model 


log g 


T e ff 


[M/H] 


log X C o 




(cgs) 


(K) 






A 


4.5 


870 


-0.5 


-4.0 ±0.25 


B 


5.0 


940 


-0.3 


-3.8 ±0.3 


C 


5.5 


1030 


-0.1 


-3.5 ±0.3 
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Table 4 

Near Infrared Colors of GL 229B 



Model logy T cS [M/H] J-K H-K K - V 



A 


4.5 870 


-0.5 


-0.11 


-0.27 


2.25 


B 


5.0 940 


-0.3 


-0.14 


-0.30 


2.26 


C 


5.5 1030 


-0.1 


-0.08 


-0.24 


2.17 




Leggett et al. 1999 




-0.10 ±0.07 


-0.07 ±0.07 


2.24 ±0.11 
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2 4 6 8 

log P (dyn/cm 2 ) 




0.9 0.95 1 

A (fim) 



0.8 - 



$ 0.6 h 



Si 

go 0.4 h 



0.2 - 




Gl 229B 
[M/H] = 
[M/H] = -0.4 




J I I I I I I I I I I I I I L 



n 1 1 1 1 1 1 1 1 1 1 1 1 1 r 



n 1 r 



0.8 - 



£ 0.6 h 



go 0.4 - 



0.2 




Gl 229B 
[M/H] = -0.2 



J I I I I I I I I I I I I I L 



J I L 



1.2 



1.22 



1.24 1.26 
A (/im) 



1.28 



1.3 



3 




chemical equilibrium CO 

all C into CO 

j I i i i I i i i 

1.56 1.58 1.6 

A (/mi) 




2.05 2.1 2.15 

A (yttm) 



[M/H] = 
[M/H] = -0.5 



Gl 229B 
[M/H] = -0.3 



2.06 2.08 2.1 

A (/im) 



i 1 1 1 1 1 1 1 1 r 




2.02 2.04 2.06 2.08 

A (/im) 




2.105 2.11 2.115 
A (yu.m) 



log P (dyn/cm 2 ) 
3 4 5 6 7 8 




log T (K) 



log P (dyn/cm 2 ) 
3 4 5 6 7 8 




log T (K) 



log P (dyn/cm 2 ) 
3 4 5 6 7 8 




log T (K) 



